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In study of pseudo(quasi)-hermitian operators, the key role is played by the positive- 
definite metric operator. It enables physical interpretation of the considered systems. 
In the article, we study the pseudo-hermitian systems with constant number of par- 
ticles in equilibrium. We show that the explicit knowledge of the metric operator is 
not essential for study of thermodynamic properties of the system. We introduce a 
. simple example where the physically relevant quantities are derived without explicit 

^ ' calculation of either metric operator or spectrum of the Hamiltonian. 
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1 Alternative representation of hermitian operators 

, In quantum mechanics, the central role is played by unitary transformations. They leave 

scalar product and norm of states unchanged. Particularly, Hamiltonians are required to 
be hermitian in order to keep unitarity of time evolution. In this section, we want to 
examine how the theory should be modified to keep the norm conserved when a wider class 
"ti ' of transformations is considered. 

1.1 Non-trivial scalar products 

j> , Let us consider a hermitian Hamiltonian H with discrete spectrum which acts on TL = 

(V, (.,.))• Let the operator B be a non-unitary bijection of TL to itself (and bijection on V 



particularly) 



gf/r 1 . (i) 

Let us use B as a similarity transformation 

H = B~ l HB, i> = B- l tj) eH, ^eH. (2) 

It leaves the vector space V unchanged. However, the Hamiltonian H ceases to be hermitian 
on TL. Instead, it satisfies the following relation 

#t = eHG~\ (3) 

where the operator is hermitian, strictly positive and has of the following explicit form 

e = b^b, e = e f , e > o. (4) 
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The norm of states is not preserved under the similarity transformation (|2|). There holds 
the following relation instead 

IHI 2 = (^,V) = W,^) = (^,e^). (5) 

To keep the norm conserved, we define an alternative scalar product (., .)e = ( . , G . ) on 
the vector space V. Consequently, the operator B may be considered as a unitary mapping 
between two Hilbert spaces 7i and 7i which coincide in vector space V but Hilbert space 
structure of Tt is induced by (., .)e- Particularly, H is hermitian in H. We may rewrite ([5]) 
as 

IM| 2 = IMIe- (6) 

To avoid possible misunderstanding caused by the alternative definition of scalar prod- 
uct, let us call the representation (H,7i) as a hermitian picture while we will call the rep- 
resentation (H,7i.) as a quasi-hermitian picture. The description of a physical system in 
both representations is equivalent as these are related by the unitary mapping B : TL i— > 7i. 

The advantage of quasi-hermitian picture lies in extended freedom in choice of B. It may 
be exploited to simplify the Hamiltonian and determine its spectral properties more easily. 
On the other hand, we lose an intuitive interpretation for the operators. For instance, the 
coordinate operator in the quasi-hermitian picture 

X = B~ X XB (7) 

may have quite complicated form dependently on the choice of B. 

Concluding, the effort to employ wider class of transformations led to the construction 
a new Hilbert space with non-trivial metric operator. In this context, we may refer to 
many-body physics where the alternative representation is used in boson realization of 
fermionic systems [HE]- The non-Hermiticity of resulting bosonic Hamiltonian is treated 
by redefinition of the scalar product which fits to the strategy outlined above. 



1.2 Statistical description of equilibrium systems 

Let H describes an ensemble of a huge but fixed number of independent particles. Let us 
focus to the statistical description of the system. The state of the system is characterized 
by a density matrix p which can be written in energy representation as 

p = ^2w n \n)(n\, Trp<oo. (8) 

n 

The coefficients W n represents (up to normalization) the probability with which we may 
find the system in the state \n). In what follows, we will consider a system in equilibrium. 
In that case, the density matrix operator is solution of Bloch equation 

%--»> < 9 » 

with initial condition p(0) = 1. Its formal solution p = e~@ H fixes the coefficients to 
be W n = exp(—/3E n ). Normalization factor of the density matrix depends on the inverse 
temperature (3 and is called partition function 

Z = Trp. (10) 
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It plays a crucial role in thermodynamics of the system as it allows direct computation of 
thermodynamic quantities. 

Let us observe how is the description modified in quasi-hermitian picture. Once we 
have © satisfying Q, the density operator in quasi-hermitian picture is 

p = B- l P B = £exp(-/3E n )|n>«ri|, (11) 

n 

where 

((n\H = E n ((h\, H\n) = E n \h), Q\h) = \n)), (h\fh)e = ((n\m) = 5 m , n 

and satisfies dpp = —Hp as long as B does not depend on (5. Let us derive the partition 
function Z for f) 1 1 1) . Partition function (1101) may be written as 

Z = ^2(n\p\n) = ^(nlBpB-^n) = ^{n\&Bp\n) = ^(n|p|n) = Z. (12) 

n n n n 

Thus the partition functions in both pictures coincide and do not depend on the mapping 
B. A similar conclusion holds for mean value of observable A; 

(Ai) P = Tr (A iP ) = Tr (i^) = (i<)p (13) 

as long as A^p is a trace-class operator. 

Bloch equation ([9]) resembles Schrodinger equation for evolution operator. The link 
can be established by the substitution (5 — > it. Then (working in quasi-hermitian picture) 
the equation ([9]) reads 

dU 

it— = HU, U(0) = l (14) 
where we defined U(t) = p(it). In x-representation, we have 

p(xi,X2)\p-nt = U(xi,x 2 ;t) = (xi\e~ lHt \x 2 ) =■ G(xi,x 2 ;t). (15) 

We should be careful with the interpretation of G. Although it is manifestly similar to a 
quantum mechanical propagator, there is a difference we should keep in mind. The vectors 
\x) are not eigenstates of the position operator in the quasi-hermitian picture. This role is 
played by \x) = B~ l \x), see (|7|). 

Despite the lack of intuitive interpretation, the "propagator" (|15p is of importance 
as it allows to compute partition function {Z = J G(x,x; —i(3)dx) with use of standard 
techniques. Particularly, let us mention standard perturbation theory or path integrals in 
this context. Dependently on the technique used, the spectrum of H is not essential for 
computation of the partition function. 

2 Pseudo-hermitian operators 
2.1 A bit of history 

Towards the end of the last century an interest of a few groups of mathematicians has 
been attracted by the rather unexpected behavior of certain manifestly non-hermitian 
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anharmonic-oscillator Hamiltonians which seemed to generate the purely real spectra O 
0]. In the next stage of development, it was Bender and Boettcher [5J [BJ who revealed 
the possible deeper relevance of the similar models (well exemplified by the "zero mass" 
Hamiltonian 

H = p 2 + ix 3 (16) 

without free parameters) in the context of physics. Subsequently, the reality of the spectra 
of the similar models has been tested [7] and, finally, rigorously proved [8]. 

As long as one has [H, PT] = in eq. (|16p . (where the operators P and T are chosen 
as representing the standard space reflection and time reversal, respectively), the authors 
of ref. [5j [6] conjectured that the reality of the spectra of the manifestly non-hermitian 
Hamiltonians of the form (]16p may be attributed to their peculiar anti-linear symmetry 
called PT-symmetry (re- interpreted as P-pseudohermiticity in ref. [9j HO]). This idea 
initiated an intensive study of pseudo-hermitian operators with anti-linear symmetries 
[111 fl~2j [131 HI"} [15] where the original specific choice of the parity is being extended to any 
invertible hermitian mapping r\ such that the operator equation 

fl1*77 = r]H (17) 

is satisfied. In such a case the spectrum of the corresponding pseudo-hermitian, diagonal- 
izable operator is either purely real or it contains the complex conjugate pairs [5[ l6l |9| [TO]. 

The operator r) is non-unique and represents rather a class of generally indefinite op- 
erators satisfying (jlT[) . However, it revealed that when the spectrum of the Hamiltonian 
is purely real, the class contains positive definite operators and, consequently, a consistent 
probabilistic interpretation (including unitary time evolution) of state vectors associated 
with pseudo-hermitian Hamiltonian is enabled [TBJ [TTJ, [TO] . With a special denotation of 
these operators, the equation (fT7|) reads now 

H*e = Off, G > 0. (18) 

Due to its positive definiteness, the metric operator can be decomposed Q = B^B. 
In this point, we may feel inspired by the previous section and say that the operators 
H and G establish quasi-hermitian picture of a physical system. The hermitian picture 
H = BHB^ 1 may provide here an intuitive insight in the nature of interaction which 
ceases to be transparent in the quasi-hermitian picture 0. Moreover, additional physical 
observables in quasi-hermitian picture may be defined with help of appropriate operators 
in hermitian picture, see eq. (J7|). 

Let us emphasize, that the current treatment differs from the previous section. The 
relation between hermitian and quasi-hermitian representation is more complicated as the 
operator Q (and B) is not given explicitly in the present case. Its construction depends on 
the Hamiltonian H and the resulting metric operator is generally non-unique. There exists 
rather a class of Q's which satisfy (|18p for given H. The construction and the associated 
ambiguity of the scalar product (., .)© represent the main questions within the theory and 
are tackled in various ways, see e.g. [U [TOJ l20| I2T]. 

The conclusions of the previous section find their important application here. We 
showed (|12p that the partition function is the same in both pictures. Thus we may compute 

1 ln (|16l) . x and p do not represent position and momentum operators. These are x = B~ xB and 
p = B~ pB in quasi-hermitian picture. 
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Z in quasi-hermitian picture where the Hamiltonian is given explicitly. This enables direct 
computation of thermodynamic properties of a physical system in equilibrium described 
by quasi-hermitian without explicit knowledge of the metric operator 0. 

In the next section, we intend to present a simple example of quasi-hermitian system. 
We derive its thermodynamic characteristics independently on both the metric operator 
and the spectrum of the Hamiltonian. 



2.2 Thermodynamics of 2 x 2 toy model 

Let us consider an ensemble of N distinguishable particles neglecting their mutual inter- 
action. The partition function may be factorized 

Z = (Z s f, (19) 

where Z s is the partition function of the subsystem. In the quasi-hermitian picture, the 
subsystem is described by the following Hamiltonian 

H = ( ^ M , a,b,eeR. (20) 

It is pseudo-hermitian with respect to r) 

H^ = r,H V -\ V=(l _°! ) (21) 

and keeps real spectrum as long as |e| < (a — b)/2. In this case, the mapping B can be 
constructed and reads 
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B= * ^ + {a z b)ib~a+d ) \ d= ^(a - 6)2 - 4e 2 . 

\ T -4e 2 +(a-6+d) 2 / 

Let us further note that the energy levels of the system may depend on the external 
parameters, e.g. on the volume v in which the particle is confined. 

We want to find Z s . It would be possible to compute the partition function exactly. 
However, let us pretend that the spectrum of the operator is hard to be found. We shall 
solve the equation (I14p . Separating the Hamiltonian into its hermitian part and non- 
hermitian perturbation term 

H = H + eV, (22) 
we introduce the following transformation 

U (t) = exp (-iH t) (23) 

and define a new operator 

Ui{t) = ul(t)U(t)U o {0). (24) 

In the interaction representation, the evolution equation (|14p of the system is given by the 
following 

i^ = V I (t)U I (t), U I (0) = 1 (25) 
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Figure 1: Graph of function A = Log 

Exp (-PEk) + Exp (-/3E 2 ) , E j = - 
improves its accuracy with decreasing e. We set a = 1, 6 = 2. 



for e = 0.03, 0.02, 0.01 where Z e s x 



+ b + (— l) 3 \J (a — b) 2 — 4e 2 ) . The approximation 



where Vj(t) = UQ(t)VUo(t). Integrating both parts of the equation, we can find the solution 
iteratively, 

U^t) = 1 + (-i)e / dhViih) - e 2 / dh / dt 2 V I {t 1 )V I {t 2 ) + 0(e 3 ). (26) 
Neglecting the higher order terms, the explicit result is 



Uj{t) 



1 a— b V 1 a— o 



1 - -S- ( it 

a—b 1 



a— 6 

-^_ e — i(a — b)t 

a—b 



(27) 



The approximative value of the partition function Z is 

Z,(J3) = Tr(Uo(t)Ui(t))\ Mf} 



„ e -*P + e -W + ±0[ e -<'< -<->>> 
a — b 



(28) 



It would be in order to discuss the range of parameters (inverse temperatur^l /3 = 1/T, 
coupling constant e) for which the approximation f|28|) is relevant. Instead, let us compare 
(I28|) with exactly computed partition function Z e s x in Fig.l. 

Now, we can compute thermodynamic quantities. The free energy of the system is 



1. 
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F = -N%LogZ s = -TV-Log ( e"^ + e~ po + e z p 

p p V a—b 



g — Pa _ e — /3b 



(29) 



and enables a direct computation of the entropy 

dF 



S 



dp- 



iVLog 



e -a(3 + e -bf3 + 



I f3 ( e -a/3 _ g-fe/3) 



"We use atomic units, i.e. k = 1. 
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[6 2 + 6(0-6 - e 2 /3)] + [-e 2 + a(g - b + e 2 /?)] 
P e a / 3 (a-6-e 2 /3) + e^(a-6 + e 2 /?) 1 j 

The internal energy of the system is 

r/ g WF) e? b [a {a - b + e 2 0) - e 2 ] + [6 (a - b - e 2 (3) + 6 2 ] 

5/3 (a - 6 - e 2 /3) + (a - b + e 2 /3) ' 1 j 

The specific heat per particle c v describes how the temperature changes when the heat is 
absorbed and the volume v of the system remains unchanged. For our system, it reads 
explicitly 

1 / dU \ _ 2 e^ a+ ^ [{a - b) 2 [(a - b) 2 - 4e 2 - e 4 f3 2 ] + 2e 4 ] - e 4 {e 2 ^ + e 2 ^) 

° v ~ N \dp ZT ) v ~ [eP a (a-b- e 2 0) + eP h (a - b + e 2 (3)f ' 

(32) 

see Fig 2. 

The pressure P of the system is defined as a derivative of free energy with respect to 
the volume v 

dv ) 1//3 



P 

9 ^e^ 2 (m (v 2 (n - m) - (3e 2 ) + v 2 e 2 ) + e^ m ' v2 (n(v 2 (n - m) + (3e 2 )-v 2 e 2 ) 



v 3 [ef 3 ^ 2 (v 2 (n - m) - (3e 2 ) + e^ m / v2 (v 2 (n - m) + j3e 2 )] 

(33) 

We assumed dependence of the energy levels in the following form 

n m e 

d — ~, ~ni £ ~n ■ 

yZ yZ yZ 

The equation ([33]) relates pressure, volume and temperature of the system and represents 
the equation of state for our system. 
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3 Discussion 

An ambition of pseudo-hermitian quantum mechanics is to serve as a consistent framework 
for description of physical processes, comparable with or exceeding the standard theory in 
predictive power. Currently, the contribution of pseudo(quasi)-hermitian operators is seen 
in the alternative representation of hermitian systems. Particularly, complicated hermitian 
Hamiltonians can be defined implicitly with help of quasi-hermitian operators which are 
easier to deal with [22], [23] . 

The construction of a mapping which relates these two operators is a crucial task of 
the theory and attracts a lot of attention pU [25j [26l E7J [28]. It is essential for proper 
probabilistic interpretation of states of non-hermitian Hamiltonian and for construction of 
additional operators representing physical observables. 

In the article, we considered statistical properties of an ensemble of pseudo-hermitian 
systems in equilibrium. Particularly, we were interested in thermodynamics of the system. 
It was the main observation that the explicit knowledge of the positive-definite metric op- 
erator G is not essential for derivation of the physical quantities. This opens a way for di- 
rect derivation of genuine physical properties of various pseudo-hermitian (PT-symmetric) 
models without tedious computation of the metric operator O (or, alternatively, Q or C 
operator, see [16} I17j. Moreover, various techniques of standard quantum theory may be 
used without any modification. In this context, let us mention path integrals [29] or opera- 
tor method [30] which might be particularly suitable for computation of partition functions 
of PT-symmetric Hamiltonians with polynomial potentials, e.g. (|16p . 
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